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Ratio tests for change point detection 

Lajos Horvath*^, Zsuzsanna Horvath^ and Marie Huskova^-^ 

University of Utah and Charles University 

Abstract: We propose new tests to detect a change in the mean of a time 
series. Like many existing tests, the new ones are based on the CUSUM process. 
Existing CUSUM tests require an estimator of a scale parameter to make them 
asymptotically distribution free under the no change null hypothesis. Even if 
the observations are independent, the estimation of the scale parameter is not 
simple since the estimator for the scale parameter should be at least consistent 
under the null as well as under the alternative. The situation is much more 
complicated in case of dependent data, where the empirical spectral density 
at is used to scale the CUSUM process. To circumvent these difficulties, new 
tests are proposed which are ratios of CUSUM functionals. We demonstrate 
the applicability of our method to detect a change in the mean when the errors 
are AR(1) and GARCH(1,1) sequences. 



1. Introduction 

Change point detection is an important part of statistical and economic analysis. 
Predictions and statistical inference will be invalid if changes in the regimes during 
the data collection period are not taken into account. The main problems in the 
change point analysis are to decide whether the statistical model for a series of 
observations does not change (no change situation) or whether the model changes 
one or more times and in the latter case to identify when the changes have occured. 
For surveys on change point methods we refer to Csorgo and Horvath [.'!] and Perron 

[11]- 

In this paper we consider at most one change in the location model 

Xk = Hk + (k i <k <n, 

where ^i, . . . , ^„ are the means of the respective observations while ei, . . . , e„ are 
random error terms with zero mean satisfying some additional assumptions specified 
below. Under the no change null hypothesis 

Ho : ^fe = /x l<A:<n 

while under the alternative 

Ha '■ there is 1 < fc* < n such that /ii = /i2 = • ■ ■ = Mfc 7^ Mfe+i = ■ • ■ = Mn- 
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The most popular methods are based on functionals of properly standardized cumu- 



lative sums (CUSUM) Ei=i(^^-^")' ^ 



I, where X„ = (1/?^) Ei<,;<„ 



For example, we reject Hq if 

k 

r„,i = max^ I ^(X, - X„)|/(na2)V2 

is large, where and ncr^ is a proper estimator for the variance of X]i<j<n Simi- 
larly, the R/S statistic proposed by Lo [10] is 



Tn.2 — 



1 



Tn.3. — 



/ 

V 



max y^(Xj - Xn) - mm y^(Xj - X„ 

l<k<n \<k<n 
. - - 1=1 - - i=\ 

The statistic r„.2 was modified by Giraitis et al. [^i] who introduced 



1 



E 

k=\ 



n A; 

fc=i i=i 



Under suitable assumptions on the error terms all three test statistics arc sensitive 
w.r.t. changc(s) in the mean (location). Asymptotic properties of Tn^i,* — 1,2,3 
were derived under the conditions 



(1.1) 

and 

(1.2) 



Eti ==0, 1 < i < oo 



there is ct > such that n^^/^ ^ e.^ aW{t) 



\<i<nt 



where {W^(t),0 < t < oo} is a Wiener process. ( ^—i' denotes weak convergence 
in X'[0, 1].) Condition (1.2) means that {e^} is a weakly dependent sequence which 
satisfies the functional central limit theorem. Since Tn,?:, i = 1, 2, 3 are functions of 
^j^^^^^j ei, (1.2) will yield the asymptotic distributions of the test statistics 

both under Hq and Ha- However, the estimator cr^ must satisfy that ct^ cr^ 
under Hq and at least it must be bounded in probability under the alternative. If 
{ei} is a strictly stationary sequence with < Ee"^ < oo, then the estimation of 
is based on the fact that it is related to the spectral density at 0. So one needs to 
choose a kernel and the number of lags used in the estimation. One of the most 
popular choices is Bartlett's estimator. However the rate of convergence is very slow 
even under Hq and cr^ might go to infinity under Ha- Modifications of the Bartlett 
estimator in the change point context can be found in Berkes et al. [1] and Berkes 
et al. [2]. Therefore it is desirable to develop procedures for testing Hq against Ha, 
where the estimator of from (1.2) is not needed. We develop such test procedures 
based on functionals of CUSUMs. 

In the definitions of Tn,i,i = 1,2,3, functionals of CUSUMs are computed for 
the first k and the last n ^ k observations. If the difference between functionals is 
large for at least one k, the null hypothesis of no change is rejected. We suggest 
computing the ratio of the CUSUM functionals instead of the differences. This way 
there will be no need for the estimation of . Instead of using T„ i we suggest 



Vn,i — max 

' n5<k<n—n5 



maxi<i<fe 






maxfc<,;<„ 
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where < (5 < 1/2 and 



i<i<fc 



and 



n ~ k 



k<i<n 



Similarly, 



Ki,2 = max 

' nS<k<n—nS 



and 



14,3 = max 

n6<k<7i—5n 



maxi<i<fc X]i<j<i(^j — -'^fe) — niini<i<fc X]i<j<i("'^i ~ -^fe) 
maxfe<,<„ Ej<j<„(^j -^fe) - minfe<,<„ Ei<j<„(^j - ^fe' 



1 

n — k 



T,'i=k+i YJj^ii^j ~ Xk) 



Om- first result gives the convergence in distribution results for Vn.i, « = 1, 2, 3 under 
Hq. Let W{t), < t < cx) be a Wiener process and define the following processes: 





= sup \W{s) - 
n<s<t 


'{s/t)W{t)\, 






= sup \W*{s) 

t<s<l 




t)W*{t)\, 


m.iit) 


~ sup {W{s) - 
a<s<t 


is/t)Wit)) - 


inf {W(s 

0<s<t 




= sup {W*{s) 




-t))W*it)) 



and 



t<s<l 

M iW*{s)-i{l^s)/{l~t))W*{t)), 

t<S<l 



VsAt) = / (Wis) - {s/t)W{t)Yds --{I (Wis) - {s/t)Wit))ds 



mAt)^ I iw*is)-i{i-s)/ii-t))w*it)rds 
1 



1 ^2 
{W*{s)-((l-,)/{l-t))W*{t))ds 



1 - 1 

where W*{t) = - 

Theorem 1.1. If Hq, (1.1) and (1.2) hold then 

mAt) 



(1.3) 



K,i sup , 

5<t<l-5 ??l,2(,Ej 



(1.4) 

and 
(1.5) 
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We note that sup5<«i_5 can be replaced with suPq^j^^^^ in (1.3)-(1.5). Since 

a.s., 



hm sup 

t-*l- t<s<l 



W*{s)-^^W*{t) 



(1.6) hm sup \ W*{s)-- — -W*{t) 

t<s<l V I -t 

- ^mf ^ [W*{s) - |^VK*(0 ) = a.s.. 



and 



1 / / l^s 



W*{s)-- W*{t)]ds] =0 a.s 



2 



we cannot replace sup^<j<]^_^- with sup5<(^]^, in (1.3)-(1.5). 

The Wiener process W has independent increments and therefore for any < i < 
1 we have that {W{s)-{s/t)W{t),Q <s<t} and {M^*(s) - ((1 - .s)/(l - t))W^*(t), 
t < s < 1} are independent. Change of variable and the scale transformation of W 
give that 

m,i{t)= sup \W{ut) -uW{t)\=t^'^ sup \B{u)\, foraU < f < 1, 

0<M<1 0<U<1 

where B{u) ~ W{u) — uW{l) is a Brownian bridge. Therefore for any < t < 1 
Vi,i{t) -D f t supo<„<i 



m,2it) \l-tj supo<„<i |B2(ii)r 

where {i?i(u),0 < it < 1} and {i?2(u),0 < w < 1} are independent Brownian 
bridges. Similar arguments give 

m,iit) V ( t supq<„<i - info<„<i -Bi(m) 



and 



mA^) yl^V supo<„<i S2(u) - info<u<i -B2(w) 
?73,i(i) 5 t !l Bl{u)du - (^jl Bi{u)du^ 



Bl{u)du - i B2{u)du 



m,2it) 1 - t 
for any <t <1. 

Kim [7] used ratio tests to detect changes in the persistence of a linear time series. 
The asymptotic as well as the finite sample properties (including power) of Kim's 
test were investigated by Kim et al. [n] and Leybourne and Taylor [0] . Since we try 
to detect changes in the means (location) our tests are different from Kim's so we 
need to investigate the asymptotic power. Let A„ = — /ifc*+i be the size of the 
change. 

Theorem 1.2. If (1.1) and (1.2) and Ha hold, 

(1.7) k* = [n9] with some < < 1, 
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and 
(1.8) 

then 

assuming that 
(1.9) 



oo, 



oo, i = 1, 2, 3, 



S <0 <1-S. 



Our tests were developed to check if the mean has changed at an unknown time. 
Ratio type tests can also be used to see if the sequence changes from "stationary" 
into "difference stationary" . We say that the sequence is "stationary" if the sum of 
the Xk's satisfies the functional central limit theorem and "difference stationary" 
if the Xk's themselves satisfy the functional central limit theorem with suitable 
normalization. Now we consider the following alternative: 



Hii : there is 1 < k* < n such that Xk 



1^1 + Ek, l<k<k*, 

^ + ffc. + . . . + efc, k*<k<n. 



The statistics Vn^i, i — 1, 2, 3 may not be able to detect if Hq or H'^ hold, since even 
under H\ the statistics have nondegenerate limit distributions. Hence we suggest 
the following modification of the test statistics to detect H^: 



ZnA 



max 

n5<k<n—n5 



niaxfc<j<„ 




-Xk) 


maxi<i<fe 




~Xk) 



z.„ 



and 



maXfc<,<„ Y.^<J<rX^3 - ^k) - minfc<,<„ 'Er<j<ni^J ' ^k) 
ITT&j^C — — - , — . — — - , — . 

nS<k<n-nS maxi<j<fe Z^i< ^^(^j - Xk) - mmi<i<k l^i<j<i(Xj - Xk) 



Zn,3 = max 

' n5<k<n—Sn 



E 



i=k+l 



EU(x,-Xk) 



1 

71— k 



Y.]=li^] - ^k) 



The limit distributions of Zn.i,i = 1,2,3 can be easily derived following the proof 
of Theorem 1.2 and we get the following results: 

Zn,i sup -— , Z„.2 sup and Z„,3 sup --. 

a<i<l-(5 '7l,U'^J S<t<l-S V2,l[t) S<t<l-S V3,l\'^) 

Theorem 1.3. If (1.1), (1.2), (1.7)-(1.9) and hold, then 

p 

Zn.i OO, i = 1,2, 3. 
1,2,3 have power against the alternative 

( H + f-k* + ■ ■ ■ + ek, l<k< k*, 



However, Vn i,i 



there is 1 < fc* < n such that X^ 



fi + Ck, k* < k < n. 
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Theorem 1.4. If (1.1), (1.2), (1.7)-(1.9) and H*^ hold, then 

p 

Vn,i ^ OO, I = 1, 2, 3. 

The alternative ff^* is somewhat the opposite of i/^; the first observations follow 
a random walk and at k* they turn into a "stationary" sequence. Of course, the 
statistics may not detect the difference between Hq and H\. If we are interested 
only if a change occured from or into a random walk at an unknown time, i/e.,we 
are testing Hq against H\ U H'^ , we must combine Vn^i and Zn^i. Let 

T„^i = max(F„^j, Z„^i) i = l,2, 3. 

Following the proof of Theorem 1.1 one can easily verify that under Ho 

Tn^i^m&yii sup j--, sup -— ^ , 

L5<t<l-5 ?7i,i(*) 5<t<i-5 V1.2V') J 

7n,2^max<^ sup --, sup ^ , 

L5<t<l-5 m,lW S<t<l-& V2,2[t) J 

Tri,3^max<^ sup j--, sup tttj ^ ■ 

l5<t<l-S V3,l[t) S<t<l-& m,2{t) ) 

Theorem 1.5. If (1.1), (1.2), (1.7)-(1.9) and hold, then 

~ p 

Tn,i CO, i = 1,2, 3. 

Remark 1.1. It would be more natural to use (5 = in our results. However, as we 
pointed out after Theorem 1.1, (1.6) yields that supg^^^j^ 7^i_i(t)/r/i,2(i) = 00 a.s. 
By Chung's law (cf. Csorgo and Rcvcsz [4]) for any v > 1/2 

lirn (1 — t)^'' 771, 2(0 = 00 a.s. 

and therefore 

(1.10) P{ sup (1 - i)"??l.l(t)/m,2(t) < 00} = 1. 

0<t<l 

By (1.10) we conjecture that 



max 1 

l<k<n \ 



Using weight functions, the statistics Vn,2 and l/„,3 can be modified in a similar 
way so one can take (5 = in the weighted statistics. 

Remark 1.2. We would like to note that ratio tests can be derived not only for 
partial sums with a Wiener limit but for more general processes. 



max I iX,Xk)\ o<t<i mAt) 
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2. Applications 

In our first example the error terms are linear processes defined as 

oo 

(2.1) efc = y^^ajSk-i, 

i=0 

where 
(2.2) 

6i, —oo < i < oo arc independent identically distributed random variables 

and 

E6o = and E6l < oo. 

If 

(2.3) |ai|<ooand ^ m ^ 0, 

0<i<oo 0<i<oo 

then (1.2) holds with cr^ = E^a{J2a<i<oo "*)^- "^^^ proof of (2.3) is in Hannan [6] 
(cf. also Wang et al. [12]). We would like to note that by Wu and Min [l::;], (1.2) 
holds for sums of linear processes without assuming (2.2). 

Wc studied the behaviour of Vn,i when 5 — 0.2. We used Monte Carlo simula- 
tions to get critical values for sup 2<t< s Vi,i{t)/'ni.2{t)- In our simulation study we 
assumed that (5^, —oo < i < oo are independent standard normal random variables 
and Ci ~ p*. This means that the observations are elements of of a stationary AR(1) 
process with parameter p. The results in Table 1 suggest that the asymptotic criti- 
cal values are acceptable even for moderate sample sizes, if p is not close to 1. The 
power in Tables 2-7 is a decreasing function of p as p tends to 1. The location of 
the time of the change has little effect on the power; the power is nearly the same 
for k* = n/2 and k* = n/4. 

In the second example we assume that are elements of a GARCH(1,1) sequence. 
This means that satisfies the recursion 

e-k — 5kTk and = + ae|_;^ -I- /3r^_;^, -oo < fc < oo, 

where a;>0,a>0,/3>0. Assuming that (2.2) holds and 

E5l < oo and aES^ + f3 < 1, 

then Berkes et al. [2] proved that (1.1) is satisfied with cr^ = uj/{l ~ aE6^ — (3). 

Table 1 

Simulated significance levels for Vn^i when S = 0.2 



P 



n 


Level 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


200 


0.1 


0.148 


0.154 


0.164 


0.183 


0.204 


0.233 


0.273 


0.338 


0.230 


200 


0.05 


0.088 


0.089 


0.108 


0.120 


0.136 


0.158 


0.204 


0.265 


0.152 


200 


0.01 


0.023 


0.024 


0.032 


0.040 


0.051 


0.069 


0.094 


0.143 


0.072 


500 


0.1 


0.116 


0.121 


0.125 


0.154 


0.152 


0.166 


0.205 


0.235 


0.326 


500 


0.05 


0.061 


0.066 


0.067 


0.078 


0.087 


0.102 


0.133 


0.157 


0.244 


500 


0.01 


0.013 


0.016 


0.017 


0.019 


0.024 


0.036 


0.051 


0.058 


0.012 


1000 


0.1 


0.120 


0.100 


0.110 


0.120 


0.130 


0.150 


0.168 


0.200 


0.230 


1000 


0.05 


0.072 


0.048 


0.052 


0.058 


0.062 


0.070 


0.096 


0.116 


0.152 


1000 


0.01 


0.016 


0.004 


0.004 


0.004 


0.008 


0.010 


0.018 


0.022 


0.072 
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Table 2 

Power of Vn.i when S = 0.2, A = 0.5 and k* = n/2 



P 



n 


Level 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


200 


0.1 


0.586 


0.528 


0.474 


0.420 


0.379 


0.366 


0.357 


0.385 


0.488 


200 


0.05 


0.467 


0.414 


0.367 


0.325 


0.289 


0.271 


0.270 


0.301 


0.411 


200 


0.01 


0.261 


0.223 


0.192 


0.165 


0.147 


0.131 


0.137 


0.166 


0.260 


500 


0.1 


0.868 


0.806 


0.720 


0.623 


0.525 


0.431 


0.359 


0.311 


0.349 


500 


0.05 


0.779 


0.669 


0.608 


0.509 


0.408 


0.326 


0.263 


0.223 


0.267 


500 


0.01 


0.551 


0.453 


0.369 


0.282 


0.212 


0.158 


0.122 


0.104 


0.137 


1000 


0.1 


0.974 


0.940 


0.892 


0.836 


0.700 


0.558 


0.424 


0.358 


0.308 


1000 


0.05 


0.940 


0.892 


0.838 


0.710 


0.592 


0.452 


0.306 


0.272 


0.216 


1000 


0.01 


0.818 


0.700 


0.596 


0.492 


0.362 


0.254 


0.148 


0.128 


0.118 



Table 3 

Power of Vn.i when 5 = 0.2, A = 1 and k* = n/2 



P 



n 


Level 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


200 


0.1 


0.962 


0.929 


0.880 


0.807 


0.718 


0.613 


0.523 


0.465 


0.511 


200 


0.05 


0.919 


0.872 


0.805 


0.720 


0.616 


0.516 


0.427 


0.381 


0.430 


200 


0.01 


0.775 


0.694 


0.604 


0.505 


0.4112 


0.325 


0.265 


0.230 


0.277 


500 


0.1 


0.999 


0.999 


0.993 


0.976 


0.935 


0.840 


0.686 


0.494 


0.406 


500 


0.05 


0.999 


0.993 


0.981 


0.947 


0.877 


0.756 


0.577 


0.398 


0.321 


500 


0.01 


0.980 


0.957 


0.905 


0.819 


0.702 


0.533 


0.356 


0.230 


0.182 


1000 


0.1 


1 


1 


1 


0.998 


0.992 


0.954 


0.846 


0.592 


0.350 


1000 


0.05 


1 


1 


1 


0.994 


0.974 


0.918 


0.756 


0.496 


0.290 


1000 


0.01 


1 


0.996 


0.988 


0.966 


0.890 


0.776 


0.536 


0.302 


0.144 



Table 4 







Power of 


when S 


= 0.2, A 


= 1.5 ai 


id k* 


n/2 






n 


Level 








P 












0.1 


0.2 


0.3 


0.4 


0.5 


II. H 


0.7 


0.8 


0.9 


200 


0.1 


0.999 


0.997 


0.989 


0.970 


0.925 


0.842 


0.719 


0.594 


0.553 


200 


0.05 


0.996 


0.987 


0.972 


0.937 


0.871 


0.764 


0.630 


0.501 


0.471 


200 


0.01 


0.969 


0.938 


0.889 


0.812 


0.706 


0.572 


0.440 


0.341 


0.327 


500 


0.1 


1 


1 


1 


1 


0.997 


0.977 


0.904 


0.717 


0.493 


500 


0.05 


1 


1 


1 


0.998 


0.989 


0.952 


0.838 


0.622 


0.402 


500 


0.01 


1 


0.999 


0.995 


0.982 


0.941 


0.837 


0.661 


0.410 


0.248 


1000 


0.1 


1 


1 


1 


1 


1 


0.998 


0.984 


0.860 


0.572 


1000 


0.05 


1 


1 


1 


1 


1 


0.998 


0.964 


0.788 


0.460 


1000 


0.01 


1 


1 


1 


1 


0.996 


0.970 


0.876 


0.606 


0.270 



Table 5 

Power of Vn.i when S = 0.2, A = 0.5 and k* = n/4 



P 



n 


Level 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


200 


0.1 


0.563 


0.512 


0.461 


0.415 


0.378 


0.350 


0.346 


0.387 


0.492 


200 


0.05 


0.443 


0.396 


0.349 


0.308 


0.275 


0.256 


0.256 


0.294 


0.414 


200 


0.01 


0.229 


0.198 


0.174 


0.149 


0.130 


0.121 


0.124 


0.159 


0.268 


500 


0.1 


0.836 


0.772 


0.687 


0.591 


0.499 


0.405 


0.341 


0.299 


0.347 


500 


0.05 


0.744 


0.662 


0.567 


0.474 


0.383 


0.308 


0.246 


0.216 


0.260 


500 


0.01 


0.496 


0.411 


0.331 


0.258 


0.195 


0.139 


0.108 


0.097 


0.138 


1000 


0.1 


0.964 


0.946 


0.894 


0.812 


0.698 


0.556 


0.412 


0.300 


0.268 


1000 


0.05 


0.942 


0.890 


0.814 


0.706 


0.580 


0.428 


0.310 


0.226 


0.176 


1000 


0.01 


0.798 


0.694 


0.592 


0.450 


0.328 


0.228 


0.150 


0.090 


0.082 
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Table 6 

Power of Vn,i when 5 = 0.2, A = 1 and k* = n/4 



P 



n 


Level 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


200 


0.1 


0.948 


0.912 


0.857 


0.785 


0.694 


0.591 


0.504 


0.463 


0.517 


200 


0.05 


0.901 


0.846 


0.781 


0.693 


0.589 


0.494 


0.415 


0.375 


0.436 


200 


0.01 


0.739 


0.659 


0.565 


0.474 


0.385 


0.308 


0.250 


0.228 


0.291 


500 


0.1 


0.999 


0.997 


0.987 


0.968 


0.911 


0.808 


0.657 


0.491 


0.404 


500 


0.05 


0.997 


0.988 


0.972 


0.927 


0.846 


0.721 


0.550 


0.391 


0.318 


500 


0.01 


0.968 


0.934 


0.873 


0.781 


0.653 


0.495 


0.344 


0.223 


0.182 


1000 


0.1 


1 


1 


1 


0.998 


0.994 


0.952 


0.820 


0.590 


0.368 


1000 


0.05 


1 


1 


0.998 


0.996 


0.980 


0.904 


0.730 


0.464 


0.278 


1000 


0.01 


0.998 


0.994 


0.962 


0.874 


0.726 


0.488 


0.266 


0.126 


0.082 












Table 7 
















Power ofVn,i 


when 5 = 0.2, A 


= 1.5 ani 


d k* = 


n/4 
















P 












n 


Level 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


200 


0.1 


0.999 


0.995 


0.984 


0.962 


0.909 


0.830 


0.712 


0.598 


0.566 


200 


0.05 


0.993 


0.983 


0.968 


0.922 


0.854 


0.753 


0.624 


0.505 


0.483 


200 


0.01 


0.961 


0.923 


0.871 


0.795 


0.682 


0.556 


0.429 


0.340 


0.343 


500 


0.1 


1 


1 


1 


0.999 


0.994 


0.966 


0.874 


0.695 


0.499 


500 


0.05 


1 


1 


0.999 


0.995 


0.982 


0.931 


0.808 


0.601 


0.402 


500 


0.01 


0.999 


0.997 


0.990 


0.966 


0.915 


0.803 


0.622 


0.409 


0.239 


1000 


0.1 


1 


1 


1 


1 


1 


0.998 


0.974 


0.826 


0.472 


200 


0.05 


1 


1 


1 


1 


1 


0.994 


0.944 


0.732 


0.386 


200 


0.01 


1 


1 


1 


1 


0.992 


0.952 


0.814 


0.522 


0.230 












Table 8 
















Power 


of Vn.i when u) = 


: I, a = 0.1 13 = 0.1 


and k* 


= n/2 






n 


Level 


A 


Sim. sign 


. lev. 


A Power 


A 


Power 


A 


Power 


200 


0.1 





0.182 




0.5 


0.615 


1 


0.968 


1.5 


0.999 


200 


0.05 





0.07 




0.5 


0.499 


1 


0.962 


1.5 


0.997 


200 


0.01 





0.02 




0.5 


0.270 


1 


0.802 


1.5 


0.975 


500 


0.1 





0.114 




0.5 


1 


1 


1 


1.5 


1 


500 


0.05 





0.06 




0.5 


1 


1 


1 


1.5 


1 


500 


0.01 





0.01 




0.5 


0.600 


1 


0.984 


1.5 


0.999 


1000 


0.1 





0.116 




0.5 


0.982 


1 


1 


1.5 


1 


1000 


0.05 





0.06 




0.5 


0.962 


1 


1 


1.5 


1 


1000 


0.01 





0.01 




0.5 


0.850 


1 


1 


1.5 


1 



The simulations arc based again on the assumption that the (5i's are independent 
standard normal variables. Comparing Tables 8 and 9, we can conclude that the 
ratio test is working well even for small sample sizes when the size of the change 
is A = (no change), A = 0.5,1,1.5. The values w = l,a = 0.1 and f3 = 0.1 
correspond to very weak dependence between the observations while the choice 
uj = 0.5, a = 0.1 and /3 = 0.7 corresponds to stronger dependence. In both cases 
the power is high and the same power was obtained for k* = n/2 and k* = n/4. 

3. Proofs 

Proof of Theorem 1.1. We can assume without loss of generality that /i = 0. Let 
Z^^i{t) = n-^/^ U and Z„,2(t) = ^ t^ 

l<.i<.nt nt<i<.n 
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Table 9 







Power 


of Vn,i, when ui = 


0.5, a = 


0.1, 13 = 


0.7 anc 


i k* = n/2 






n 


Level 


A 


Sim. sign. lev. 


A 


Power 


A 


Power 


A 


Power 


200 


0.1 





0.128 


0.5 


0.478 


1 


0.967 


1.5 


0.995 


200 


0.05 





0.071 


0.5 


0.372 


1 


0.842 


1.5 


0.985 


200 


0.01 





0.018 


0.5 


0.179 


1 


0.637 


1.5 


0.916 


500 


0.1 





0.107 


0.5 


0.775 


1 


0.998 


1.5 


1 


500 


0.05 





0.060 


0.5 


0.663 


1 


0.990 


1.5 


1 


500 


0.01 





0.012 


0.5 


0.423 


1 


0.938 


1.5 


0.998 


1000 


0.1 





0.116 


0.5 


0.954 


1 


1 


1.5 


1 


1000 


0.05 





0.064 


0.5 


0.892 


1 


1 


1.5 


1 


1000 


0.01 





0.008 


0.5 


0.744 


1 


0.998 


1.5 


1 



Condition (1.2) yields that 



(3.1) 



where W*{t) = W{1) - W{t). Since 



n sup 

l<i<nt 



l<j<i 



sup 

0<i<rit 



Zn,lii/n) - ■r—-rZn,l{t) 

[nt] 



and similarly 



n sup 



nt<i<n 

by (3.1) we have for aU < (5 < 1/2 



sup 

7it<i<n 



Znai^/n) \ Zn,2{t) 



n sup 

0<i<nt 



l<J<i 



-1/2 

, n ' sup 

nt<i<n 



z<j<n 



a{ sup \W{s)-{s/t)W{t)\, sup \W*{s)-{{l-s)/{l-t))W*{t) 

0<s<t t<s<l 



Hence the proof of (1.3) is complete. The statistics Vn,2 and Vn.3 are also continuous 
functionals of Z„_i(i), Z„_2(0i < t < 1. Hence the arguments in the proof of (1.3) 
can be repeated. □ 

Proof of Theorem 1.2. Let k > k* . Then the definition of Xj gives 



^(X„ - X,) 

i=i 



j=l '^3 k 2^J = 1 



i{k-k*) 



if 1 < i < fc* 



E-=i^.-iE-=ie. + (*-fc*)A„ 

If fc = [nr] with some 6 < t < 1 — S, wc get that 



iik-k") 



k 



if k* < k. 



n '^1'^ max 

l<i<fc 



= Op(l) 



i<i<fc* 
k*{k - k*) 
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Since there is no change in the means of Xk, Xk+i ■ ■ ■ , Xn, by Theorem 1.1 we have 
that 



-1/2 

n ' max 

k<i<n 



i'^j^n 



V 



cr(l-r)i/2 sup \B{t)\ 
0<t<l 



where < t < 1} is a Brownian bridge. Observing that k*{k-k*)n-^l'^\I\n\l 

fc ^ oo we conclude that Vn^i — > c». Similar arguments yield the proof when i = 2 
and 3. □ 

Proof of Theorem 1.3. It follows from condition (1.2) that 



n ^1"^ max 
i<i<fe* 



\<3<i 



, n ' max 

k* <i<n 



E (^.-^^ 



(T sup 



sup 

9<t<l 



W{t) - -Wi9) 



(W{s) - W{0))ds 



1 - t 



{w{s) - w{e))ds 



p . . . p 

proving that Z„^i —> oo. Similar arguments give that Zn,i — > oo when i = 2 and 3. 

□ 

Proof of Theorem I.4. Simple modifications of of the proof of Theorem 1.3 gives 
the results. □ 



Proof of Theorem 1.5. It is an immediate consequence of Theorems 1.3 and 1.4. □ 
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